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SECTION  I 


INTRODUCTION 


One  of  the  most  powerful  methods  of  measuring  frequency  of  several 
simultaneous  signals,  with  limited  frequency  resolution  required,  is  by 
using  a bank  of  contiguous  filters.  If  the  incoming  signal  is  continu- 
ous wave  (CW) , the  frequency  is  relatively  easy  to  determine,  since  the 
insertion  loss  characteristics  of  the  filter  are  known.  The  relative 
outputs  of  two  or  more  adjacent  filters  then  identify  frequency.  If  the 
incoming  signal  is  pulse  modulated,  the  output  is  a complex  function  of 
time,  and  requires  more  sophisticated  detection  schemes.  To  optimize  the 
filter  bank  and  the  detection  scheme,  one  has  to  know  the  output  from 
different  types  of  filters.  It  is  time  consuming  and  expensive  to  build 
the  circuits  and  measure  their  performance  (Reference  1). 

The  purpose  of  this  report  is  to  present  an  analytical  procedure 
with  computer  calculations  which  determines  the  time  response  of  the 
output  signals,  given  the  filter  transfer  characteristics.  With  this 
method,  a large  number  of  filters  can  be  analyzed  for  the  design  of  a 
channelized  receiver  without  going  through  the  construction  and  test 
phase.  To  document  the  validity  of  the  analytical  calculations,  they 
are  compared  with  measured  performance  of  various  filters. 

The  computer  program  is  written  to  handle  several  kinds  of  band-pass 
filters,  including  Butterworth,  Chebyshev,  and  a combination  of  two  of  these 
in  cascade.  The  three  kinds  of  pulsed  carrier  input  signals  used  in  the 
calculations  are  a square,  a composite,  and  a sine  square.  The  square 
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pulse  is  assumed  to  have  infinitely  sharp  leading  and  trailing  edges.  The 
composite  pulse  has  finite  leading  and  trailing  edges  that  can  be  selected 
independently.  The  sine  square  pulse  is  commonly  used  to  approximate  a 
short  signal  with  slow  rising  and  trailing  edges. 

SECTION  II 
MATHEMATICAL  MODEL 

In  the  calculations,  the  transfer  functions  of  the  filters  are  first 
predicted  from  the  insertion  loss  versus  frequency  measurements;  then  the 
time  domain  impulse  response  of  the  transfer  function  is  derived.  The  com- 
puter program  then  calculates  the  convolution  of  the  input  signals  and  the 
filter  transfer  function. 

1.  FILTERS 

The  filters  are  either  specified  or  measured  and  their  transfer  func- 
tion analytically  fitted  to  the  data,  using  a Butterworth  or  a Chebyshev 
model.  The  basic  equation  used  to  determine  the  number  of  poles  is 


where  V is  the  output  voltage  of  the  transfer  filter  function,  w is  the  angu- 
lar velocity,  and  ui0  is  the  center  angular  velocity  of  the  filter.  For  a But- 
terworth  filter,  the  poles  are  located  as  follows: 

PR  = Au>  [cos  tt)  + j sin  n)]  (2) 

where 

k **  1,2, . . . ,2n, 

Am  * 2w  x bandwidth  of  the  filter, 
n ■ number  of  poles  of  the  filters. 
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Only  with  positive  real  parts  are  kept,  and  the  poles  corresponding 
to  a band-pass  filter  are 
Pu 


sk  - -r  * 


(3) 


For  the  Chebyshev  filters. 


Pk  “ Aw  (tanh  v sin  uk  + jcosuk) 


uk 


(2k-l) 


v • — sinh-1  ( — ) 
n ' £ ' 


(A) 

(5) 

(6) 


and 


c - Vl°R£/1°  - 1 

where  k * l,2,...,2n, 

n » number  of  poles  of  the  filters, 
Rf  - ripple  factor  in  dB, 
w = 2n  x signal  frequency,  and 
w0  - 2 ti  x filter  center  frequency. 


(7) 


Similarly,  only  Pk  with  positive  real  part  are  kept  and  the  locations  of 
the  poles  are  calculated  by  equation  (3).  The  transfer  function  of  the 
filters  is  given  as 


H(s) 


(Aw)nSn 


(S+S1)(S+S2)...(S+S2n) 


A i A2  A2n 

+ + . . . + 


S+Si  s+s2 


S+S 


2n 


(8) 


T 


and  the  corresponding  Impulse  response  is 

H(t)  ■ Aje  + Aje  ^*>t  + ...  + A2ne  ^2nt  (9) 

2.  SIGNALS 

The  general  shape  of  the  input  signal  can  be  written  as 


R0  (t)  f i u> t -1u,tl 
R(t)  - Ro (t ) cos  u)t  * — r—  |e  + e J 

where  u ■ 2*  * the  carrier  frequency  of  the  signal  and 
of  the  input  signal.  For  a square  wave  envelope  input 


(10) 

Rq  (t ) «■  envelope 
from  t » 0 to  t “ t i , 


Rq  (t ) - U(t)  - U(t-tj) 


(ID 


where  U(t)  is  the  step  function.  For  a composite  signal  as  shown  in 
Figure  1, 


R0(t)  = DtU(t)  - D(t-ti ) U ( t — t j ) - E(t-t2)U(t-t2)  + 

E(t-t3)U(t-t3) 


where  D 


— and  E 
tl 


1 

t3-t2 


For  a sine  square  input. 


(12) 


R0(t)  - sin2  — U(t)  - sin2  -1  U(t-tx ) 
tj  ti 


(13) 
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SECTION  III 


OUTPUT 


The  output  from  the  filter  is  given  by  the  convolution  integral 


Y(t)  - R(r)H(t-T )dx 


1.  SQUARE  INPUT  SIGNAL 

For  the  square  input,  Equation  14  becomes 


ft  ft 

’(t)  - U(t)  R(x)H(t-x)dx  - U (t-t i ) R(x)H(t-x)dx  (15) 

' O 


For  t < t\ 


T [vk(t)  - YU(0>] 


and  for  t > .ti 


Y(t)  = f [Yk(t,)  - Yk(0)] 


where 


Ake-S''t  re6'''  „Ck' 
Vk<T> 2 'V+C1T 


Bk  * Sk  + Ju 


ck  “ Sk  * 


" 


2.  COMPOSITE  SIGNAL 

For  the  composite  signal,  the  output  for  t < tj  becomes 
2n 


Y°i(t)  = l [Yxk(t)  - DXk(0)J 


(19) 


where 

Xk(T> 

for  tj  < t < t2 


^ e"^  L1*1/1 — \+  eCkY± i-V 

W Bk2j  (ck  CkY 


(20) 


2n 

Y (t)  = l -DXk(0)  + DXk(tl)  + DtiYk(t)  - DtjYk(ti ) (21) 

k=l  L J 


for  t2  < t < t3 


2n 

Y(t)  = l -DXk(0)  + DXk ( t j ) + Dt]Yk(t)  - DtjYk(ti)  - 
k=l  L 

EXk(t)  + EXk(t2)  + Et2Yk(t)  - Et2Yk(t2)J  (22) 


and  for  t3  < t 


Y (t)  « l r-DXk(0)  + DXk(tl)  + Dt!Yk(t)  - Dt1Yk(t1)  + 
k=l  L 

EXR ( 1 2 ) + Et2YR(t)  - Et2Yk(t2)  - EXk ( 1 3 ) - 
Et3Yk(t)  + Et3Yk(t3)]  (23) 
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J.  SINK  SQUARE  WAVE  SIGNAL 

For  the  sine  square  wave,  the  output  for  t < tj  becomes 
2n 

Y(t)  - l 
k-1 


[vk(0  - Yk<°>] 


(24) 


and  for  t < t i 


where 


2n 

Y(t)  ■=  l 
k-1 


Ok(ti>  - zk<°>] 


(25) 


A computer  program  which  calculates  Equations  16  through  26  and  plots  the 
results  has  been  written  and  is  shown  in  the  appendix.  The  following 
section  of  this  report  discusses  the  computer  program. 
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SECTION  IV 
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COMPUTER  PROGRAM 

The  flow  diagram  of  the  computer  program  shown  in  Figure  2 is 
presently  set  up  to  handle  a maximum  of  n « 12  poles,  hut  with  minor 
changes  it  can  handle  a larger  number  of  poles.  Since  only  the  envelope 
of  the  output  Is  of  interest,  the  program  plots  the  envelope  rather  than 
the  RF  signal  inside  the  pulse.  The  program  will  pick  the  minimum  of  2 n /to 
and  2 it / ut0  and  use  it  as  a standard  time  value.  The  first  thirty  local 
■naximums  are  found  using  steps  of  l/10n  times  the  standard  value.  In  order 
to  minimize  computation  time,  subsequent  points  are  found  by  incrementing 
0.8  of  the  standard  value  from  the  last  maximum  and  then  using  steps  of 
1/I0n  times  the  standard  value  to  find  the  next  maximvim.  The  estimated 
error  in  amplitude  is 


1 


cos 


S( 


1 360  \ 

2 ’ lOn  / 


0.005 


or  0.5Z  with  respect  to  the  true  maximum  value. 


SECTION  V 

CALCULATED  AND  EXPERIMENTAL  RESULTS  FOR  CONVENTIONAL  FILTERS 
In  order  to  cover  all  the  possibilities,  seven  cases  are  calculated 
with  the  results  as  shown  in  Figure  3a  through  Figure  5.  In  these  Figures, 
only  the  envelopes  of  the  signals  are  plotted  and  their  amplitudes  are 
normalized  according  to  the  peak  value  of  each  individual  output.  In 
Figures  3a,  3b, and  3c,  a 3-pole  Chebyshev  filter  is  used  with  a center 
frequency  of  300.8  MHz,  bandwidth  of  17.8  MHz  and  a ripple  factor  of 
0.25  dB.  The  input  signals  are  all  the  composite  waveform  with  Tj  - 15  ns, 
T?  - 190  ns  and  T3  - 200  ns.  The  center  frequencies  of  the  input  used  are 
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301  MHz,  309.8  MHz,  and  317  MHz,  respectively.  In  Figures  4a,  4b,  and 
4c,  a 5-pole  Butterworth  filter  is  used  having  a center  frequency  of 
300.5  MHz  and  bandwidth  of  15.8  MHz.  The  input  signals  are  all  square 
wave  with  a pulse  width  T of  300  ns.  The  center  frequencies  of  the  input 
signals  are  300.5  MHz,  291.4  MHz,  and  292.6  MHz,  respectively.  In 
Figure  5,  the  same  5-pole  Butterworth  filter  is  used.  The  input  signal 
is  a sine  square  with  Tj  - 20  ns  and  center  frequency  of  300.5  MHz.  The 
amplitudes  of  all  the  outputs  are  normalized  to  unity  and  the  times  used 
are  either  50  ns  or  100  ns  as  shown  in  the  figures.  For  a verification 
of  the  calculated  results,  experimental  measurements  are  made  for  the 
same  input  waveforms.  The  RF  pulse  was  generated  by  a RF  switch  and  the 
output  displayed  on  a Tectronix  type  7904  oscilloscope. 

Figures  6 and  7 show  the  insertion  loss  vs  frequency  of  the  throe 
filters  used  for  the  analytical  calculations  above.  The  input  signals 
are  shown  in  Figures  8,  9,  and  10.  Figure  8 shows  the  composite  signal, 
with  Tj  = 15  ns,  T2  ■ 190  ns, and  T3  = 200  ns;  Figure  9 shows  the  square 
wave  signal  of  T * 300  ns;  Figure  10  is  the  sine  square  signal  of  20  ns. 
The  outputs  are  shown  in  Figures  11a  through  13.  The  frequencies  are  the 
same  as  specified  in  Figures  3a  through  5.  Figures  11a  through  12a  match 
the  calculated  results  very  well.  However,  the  results  of  Figures  12b  and 
12c  do  not  match  Figures  4b  and  4c  as  closely.  By  changing  the  frequency 
from  291.4  to  292.86  MHz  and  from  292.6  to  293.444  MHz  as  shown  in 
Figures  12d  and  12e,  they  match  the  calculated  results  better.  Obviously, 
the  frequency  measurement,  or  the  bandwidth  of  the  filter  measured,  is 
off  slightly.  At  the  center  of  the  filter  the  shape  of  the  output  signal 
is  not  very  frequency  dependent,  but  at  the  edge  of  the  filter  the 
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output  depends  very  strongly  on  the  frequency.  In  Figure  13,  the  measured 
result  agrees  with  the  calculated  results  rather  well.  However,  the 
Input  signal,  as  shown  In  Figure  10,  has  some  reflections  in  the  time 
domain;  and  Is  not  a true  sine  square  signal.  Thus,  the  output  may 
deviate  slightly  from  the  ideal  case. 


SECTION  VI 

APPROXIMATION  OF  SURFACE  ACOUSTIC  WAVE  (SAW)  FILTER 
A SAW  filter  with  an  insertion  loss  vs  frequency  as  shown  In  Fig.  14 
is  used  as  the  filter.  This  SAW  filter  was  built  by  a cosine  square  time 
domain  configuration  on  a pedestal  (Ref.  3).  This  filter  cannot  be 
properly  approximated  by  a single  Butterworth  or  Chebyshev  filter.  However, 
it  can  be  represented  by  two  conventional  filters  in  cascade.  Both  filters 
used  in  the  analysis  are  Butterworth  with  the  same  center  frequency  of 
343.5  MHz.  One  is  a 7-pole  filter  with  bandwidth  of  9.5  MHz,  the  other 
is  a single  pole  with  bandwidth  of  7.0  MHz.  The  input  signal  is  a composite 
one  with  Tj  - 25  ns,  T2  **  320  ns,  and  T3  - 340  ns.  The  center  frequencies 
are  343.5  MHz,  335.6  MHz,  and  353  MHz.  The  calculated  and  measured  results 
are  shown  in  Figures  15a,  b,  c,  and  16a,  b,  c,  respectively.  In  Figure  15a, 
c,  and  16a,  c,  the  results  match  fairly  well;  however,  in  Figure  15b  and 
16b  the  relative  amplitudes  between  the  two  peaks  are  reversed  in  the 
calculated  and  measured  results.  It  is  suspected  that  the  approximation 
is  not  close  enough  and  may  require  a better  combination  of  cascaded 
filters.  For  some  SAW  filters,  the  analytical  model  using  Butterworth 
or  Chebyshev  filters  must  be  modified  by  including  the  phase  response  in 
order  to  obtain  the  true  time  domain  performance.  A number  of  SAW  filters 
are  presently  being  measured  to  completely  characterize  their  frequency 
domain  performance  and  subsequently  their  time  domain  response. 
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SECTION  VII 


CONCLUSION 

For  the  limited  number  of  cases  illustrated  above,  the  computer 
program  does  generate  a very  accurate  result  for  use  in  predicting  the 
output  of  a filter.  This  program  can  be  used  to  design  receivers  with 
filters,  especially  a channelized  receiver  which  uses  many  filters. 
This  program  can  also  be  extended  to  handle  low-pass,  high-pass,  and 
band-rejection  filters  by  using  the  appropriate  filter  transfer 
characteristics  (Reference  2). 

It  also  demonstrated  that  the  transient  response  of  some  SAW 
filters  can  be  approximated  by  this  approach.  The  particular  SAW  fil- 
ters that  may  be  used  in  this  approach  are  those  displaying  linear 
phase  response  in  addition  to  providing  Butterworth  or  Chebyshev 
amplitude  response  in  the  frequency  domain. 


REFERENCES 

1.  Hollis,  R.,  "Double  Detection  Filter  Techniques  Test  Report," 
Watkins-Johnson  Co.,  Palo  Alto,  CA,  1973. 

2.  Chausi,  M.S.,  "Principles  and  Design  of  Linear  Active  Circuits," 
Chapters  4,  15,  McGraw-Hill,  Inc.,  1965. 

3.  Szabo,  T.L.,  and  Slobodnik,  A.J.,  Jr.,  "Diffraction  Compensation  and 
Periodic  Apodized  Acoustic  Surface  Wave  Filters,"  IEEE  Trans,  on  Sonics 
and  Ultrasonics  U.  SU-21,  April  1974,  pp  114-119. 
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OPTIONS 


I FILTER  BWTH 
RIPPLE  FACT .1 


| FILTER  FREQ. 
f NR.  OF  POLES 
f INPUT  TIME 
SIGNAL  FREQ.) 


CALCULATE 
"P"  VALUES 
FOR  BUTTER- 
WORTH  FILTER 


_>✓''(  TYPE  OF 
"Icrv  FILTERl"2o" 


CALCULATE 

"S" 

VALUES 


CALCULATE 
"A"  & "H" 
VALUES 


LISTING 


CALCULATE 
"P"  VALUES 
FOR  CHEBYSHEV 
FILTER 


Fig.  2 Computer  program  functional  flow  diagram 
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m 


▼ 


NORMALIZE 

ENVELOPE 


CALCULATE 

PLOTTER 

COORDINATES 


Fig.  2 Computer  program  functional  flow  diagram  (Con't) 
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300.8  MHz 


FT 

I i 


Fig.  3 Calculated  output  of  a 3-pole  Chebyshev  filter  with  fo  “ 

f - 17.8  MHz,  RF  - 0.25  dB,  the  input  signal  is  composite  with 
Tj  ■ 15  ns,  T2  ■ 190  ns  and  T3  ■ 200  ns 


o.oo  0.0 « 0.10  0.1$  o.ro  o.n  o.yo  o.»  o.  *o  o. «»  o.  so 

r 1 ne  iseconos)  «to* 


3 a.  Input  frequency  301.1  MHz 


Fig.  4 Calculated  output  of  a 5-pole  Butterworth  filter  with  fo  “ 300.5  MHz 
f ■ 15.8  MHz  and  square  input  signal  of  T * 300  ns 


s 


4 a.  Input  frequency  300.5  MHz 


I 

i 


- 


300  MHz  20  MHz/d  Iv  10  dli/div 


6 Insertion  loss  vs  frequency  of  n 3-pole  Chebyshev  filter 


f0  * 300  MHz  20  MHz/d  lv  10  dB/div 
7 Insertion  loss  vs  frequency  of  ;i  S-pole  Rutterwortb  filter 


I i 

I . A 

■J 

m ■ 

i 

1 

1 

1 

I- 

■ : m 

V- 

I 

»» 

• 

M 

WL 

J 

Fig.  10  Sine  squared  signal  T = 20  ns 


11  a.  Input  frequency  301  MHz 


22 


Fig.  15  Calculated  output  from  the  SAW  filter 


TIME  iniCfOSECONOS) 


Input  frequency  343.5  MHz 


T I PIE  miCROSECONOSl 


Input  frequency  335.6  MHz 
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APPENDIX 


COMPUTER  PROGRAM  LISTING 


PROGRAM  FILTER  ( INPUT  , OUTPUT,  TA=>E4,  T A PE5  = IN°JT  , TA  PE6=OUTOUT> 

C 

C THIS  PPOGRAM  PROCESSES  AN  INPUT  SIGNAL  THROUGH  AN  N POLE  FILTER. 

C THE  OUTPUT  SIGNAL  MAS  A RFLATIVE  AM»lITUDE  AS  A ="JNCTION  OF  TIMe 
C THe  PROGRAM  IS  WRITTEN  IN  FORTpaN  IV  EXTENOFO  c3R  a CQC6600. 

C PLOTTING  IS  DONE  ON  AN  OFFLINE  CAICOMP  PLOTTER. 

C 

COMPLEX  X1,X2,X3,*4,XB,X6,X7,X8,  X88 
COMPLEX  04 

COMPLEX  X61,X62,X63,X71,X72,X73 

COMPLEX  A (4  8)  ,0(1*8)  , C (46)  , 0ME5  AI  , OH E GA  J , ° ( 1 2)  , S(  48)  , ST  , ST  2 , T FRM  4 , 
1Y  (3)  , SS ( 24) 

DOUBLE  DFLT,T 
LOGICAL  IY,IE,IFA,IENO 
LOGICAL  IFP 

DIHF  NSIOM  AH (20 0 0)  , DATA  < 10  2 4)  ,TH  (2000 
OIMFNSION  FF(S) , Sr (5) 

CDMMON/VAPS/AH, AM°T , DELT , I V AP ,L,M,T,TH,TMAX,XVAR,XVAR1, 
1YVAP,Y,OFLT2,IFO, IF A 

NAMFLIST/VAL'JF/KIS,SCF,Tl,T2,T3t  KF,  ,N P,  FCF , 3W , RF,  TMAX  , C T , I V A R , XV  A R , 
1 YVAR,IY,IE,IFA,IENn,ICF,IFP 
NAMFLIST/FIL2/KF,-»W,N°,RF,FrF 

DATA  IVAPtIFNn,TVf IE,IFA,TMAX,RI , XV AR, YV AR/ 0 , . T. , 3* . F . , 0 . , 

1 7.14153265358970,7.75,2.5/ 

SINHI(X)  =ALOG(X«-SPRT<X**?*l  .)) 

ICF=1 


EACH  FILTER  SIMULAT'D  MUST  BE  S°EDIFIED  BY  A /VALJE/  NAMELIST. 
NAMELIST  7FIL27  CAN  ONLY  BE  SPFDI-IEO  rDR  A DOH^OSITE  FILTER  AND 
MUST  IMMEDIATELY  FOLLOW  ITS  RELATED  FILTER'S  /VALUE/  INPUT. 

THE  TNPUT  OAPAMCTFRS  ftRc  INP'JTFO  THROUGH  N AMEL IST/ V ALUE  / AND  ARF.  AS 
FOLLOWS  t 

*TS  = “KINO  OF  INPUT  SIGNAL",  MAY  HA Vc  A VALJE  FROM  1 THROUGH  3 
IF  1,  THEN  THE  SIGNAL  IS  A SQUARE  WAVE 

IF  2,  THcN  A PAMO  TO  °EAK,  PEAK  for  SOME  TIME,  AND  THEN  A 
NEGATIVE  SLOPE  TO  C AMPLITUDE  AT  SOME  LATER  TIME 
IF  3,  THEN  THE  WAVP  IS  APPROXIMATED  BY  SINE  SOUAREO 
SCF  - “SIGNAL  CENTER  FREOUSN3Y"  IN  HFRT7 

T1  = TIME  DUPATION  OF  SIGNAL  FOR  Kis  = 1 OR  3 AND  FOR  <IS  = 2 
T T IS  THF  END  0?  THE  RAM*  (PEAK) 

T2  = TIME  DURATION  OF  =>EAK  SIGNAL  f T1  for  <js  = 2 ONLY 

T*  = TIME  DURATION  OF  NEGATIVE  SLOPE  * Tl  + T2  FOR  KIS  = 2 

ONLY 

KF  r "KIND  OF  FILTFR",  MAX  HAVE  A VALUE  OF  1 OR  2 
IF  l,  THrN  THE  FILTER  IS  A pUTTERWDRTH 


IF  ?,  THCN  THE  FILTLR  IS  A THEBYSHEV 
Mp  - NUMBER  OF  POLES  FOR  THE  FT  L TER  (PRESENTLY  LIMITED  TO  121 
FCE  = “FILTER  CENTER  ERE OUENO Y"  IN  HERT7 

BW  = PANOWIOTH  OF  THE  FILTER  (3EtWFEN  -3  OB  POINTS)  I 

RF  = RIPPLE  FACTOR  FDR  CH-BfSHEV  FILTFR 

TMAX  = maximum  TIME  DURATION  OF  OUTPUT  SIGNAL  (SEE  NOTE  BELOW) 
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X1IS  PACE  IS  BEST  QUALITY  PRACTICABLE 
PROM  OOP!  FURNISHED  TO DAQ 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

c 

c 

c 

c 

r, : 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r- 


c r 


TVAR 

XVAP 


WAR  = 

icf=i 

TCP-? 

i y 


ir 


if® 

i pa 


I Fwn 


M0Tc 


(INCHES)?  DDES  NOT  CHANGE 


OUT 

Y ARRAY 


MULTIPLIER  WHICH  IS  THE  NJHBER  OP  TIMES  THE  OUTPUT 
SIGNAL  EXCEEOS  THE  INPUT  SIGNAL  (SEE  NOTE  BELOW) 
NUMBER  0^  ITERATIONS  THAT  ARE  USED  IN  CALCULATING  Y(I) 
tENGTH  nc  X-AXIS  OF  °LD  T (INCHFS)  INITIALLY?  THIS  VALUE 
IS  RECALCULATED  IP  THE  A H C I ) OP  AL(I)  ARRAYS  ARE  TOO 
LARGE 

HEIGHT  Or  Y-AXIS  OP  PL  D T 
SINGLE  filter 
COMDOSI-E  FILTFR 

= DECISION  WHETHER  TO  WRITE  OUT  Y ARRAY 
IF  .TRUE.  THEN  T Hp  ARRAY  IS  WRITTrN  OUT 
IF  .FALSE.  THEN  THE  ARRAY  IS  NOT  WRITTEN 
= PECISTON  WHETHER  TO  WRITE  OUT  LNVEL OaS  Dc 
IP  .TRUE.  THEN  THF  APRAY  IS  WRITTFN  OUT 
IF  .FALSE.  THEN  THE  ARRAY  IS  NOT  WRITTEN  OJT 

=FREDUENCY  PLOT PLOT  IS  0 ONE  IF  .TRUE. 

: DECISION  WHETHER  to  WRITE  OUT  FREQUENCY  ARRAY 
IP  .T.  THFN  ARRAY  IS  PRINTED 
IF  .F.  THFN  ARRAY  IS  NOT  PRINTED 
= OFCTSTON  WHETHER  THERE  IS  ANOTHER  SET  Or  DATA 
IE  .TRUE.  THEN  A HOT  Hf  R S LT  IS  RUN 
IF  .FALSE.  THcN  THF  PROGRAM  TERMINATES 
IF  THAX  IS  SET  AT  C»  THEN  T^AX  = AH  AX  1 ( T 1 , T 3) * CT , ANO 
IE  I MAP  IS  SET  AT  0,  THEN  IVAR  = A M AX  1 ( FCF , SC F) *2. * PI * 
P3.*TM*XH  . 

FOR  COMPOSITE  FILTERS,  PARAMETERS  FOR  THE  SECOND  fjltER  ARE  INPUT 
THROUGH 

<F 
BW 

PARAMETERS. 

FORMAT 


TO  BE  RUN 


COL 
COL 
r ol 
roL 


NAME 

LIST  /FTL27 

ANO 

A 

PE  AS 

FD 

LLOWSt 

no  FCC 

RF  SEE 

N A Mr 

L 

1ST  f 

v/AL 

UE.7  FOP 

RS. 

r NA 

MELIST  INPUT 

is: 

1 

BL  ANY 

2-7 

•IVALIJF*  OP 

*?E 

I 

L2  * 

EXPLANATION  FO 


IF  omf 
CAROS 


3 BLANK 

B-BC  PARAMETER  NA**E  FOLLOWED  nY 

parameter  VALUE.  SEdERATE 
CNn  WITH  A DOLLAR  SIGN, 

CAPO  IS  INSUFFIC IrNT  - IN°J  T MAY  BE 
(FTAPT  TN  col  2),  BREAKING  9 ARA METER 


FOUAL  SIGN 

parameters 


FOLLOWED  BY 
0 Y COMMAS  ANO 


CONTINUED  ON  OTHER 
LIST  AT  ANY  COMMA. 


CALL 

PLOTS 

(OAT 

A, 

1024, 

CALL 

PL  0T  ( 

0 . ,F. 

i * 

WRIT 

r 

(*,2 

BCD) 

TSTART 

= SFC 

ONO  (0 

°> 

o-fln 

IF 

,VAL 

IJF) 

FF(1 

) * 

KF 

*EF( 

’) 

-BW 

CONT 

I" 

'|F 

IE(N 

P. 

GT.l 

2)  GO 

T 

0 7 J 0 

WPTT 

r 

(F,  V 

ALUE) 

0*1  = 

?. 

FOP 

0 MCG 

AI 

=CMPLX (0. 

,0 

“I) 

£ Fc ( 3 ) =N9 


IFF (4) =RF  ?FF(5)=FCF 


34 


.inn  d o n r>or>  -in  r> 


THIS  PAGE  IS  BEST  QUALITY  mCTJLCMU& 
EWUii  Cm  i E,UWUSHiX>  IVDftQ  - • ' 


0“J=?.»f*T*Srr 
OMFT,  AJrCMPLVn.  ,0*«J) 
NP2  = Nr>*  ■> 
TrRHirnu»,i,*oT 

MtO 

GO  Tfi  (10,  ?o»  , *F 


1C  00  11  I =1 ,NP 

TFPM2=FLOf  T t ?*I  -1  »NP)  /FL  OS  T (Mi?)  *P[ 
T rp  = --f  RMl  *COS( T^rm?) 
TPIrTFRMi»STN(Tcc>‘1?) 

°(I)  =FM"LX  < T(JR,  THf) 

11  w PI  t «-  <«, , 3090)  I , n(I  ) 

rn  to  m 


CALCULATION  Oc  VALUES  FO?  C H i 3 Y S '•* E V F I L T r ? 


?C  E I = 1 ./SORT  (1 0.**  ( >F/1 j.» -1 . > 
V*l. /FLOAT (WO)* SI MHI f - I » 

no  21  1=1 ,Np 

u=pi/float(np?>  *float  (2*i-i) 
TPPzTAMH(V)  *SIM  (l))*TF  ? M 1 
Tni=C0S<ll)#TE9Ml 
P ( T)  =rMr,LX(TPP,TnI) 

?1  wpitc  ((,,  30  011  I,n(I) 


CAlCtll  ATI  ON  nF  VAL'IFO 


30  J =0 

"0  31  I = ?* 

J = .IH 
1 **1  = I *■  1 

c(t-d  =p(.m  / ? « ♦ oh *r,e i 
F (I) =P(  I)  /•>.  -OMrGAJ 

ti  wpitf  <s,3ori*  i mi  ,s<i-n  , i (i) 

GO  TO  ) ic  f 


CALCULATION  OF  "A"  VALUrS 


CONT  INUF 

no  32  1 = 1, ‘IP? 

ST=-sm 

TFRMUrGMPl  x ( 1.  , 0.  ) 

on  33  J=1,MP? 

IF  TI.EO.JT  GO  TO  33 
TfRMA  = TFPMA* (ST  *S I J)  ) 
33  contimuf 
ST?=ST**NP 


CALCULATION  OF  V3LUFS  FOR  OllTT  I R W3  R T M FHT-* 


Si?  14  B*ST  QUAWW  HUcriCABL* 
oapy  PUKiusajj)  rox^c  , 


r WRITF  ('>,’Cei>wST’  = -,ST2,,,TFRK,=»,  TrR-4 
2C31  FCPMaT(lX,A4,?f  19.  ir,  5X,  to,  2?  IP.  1C* 

A (II  :ST’/TC9% 

32  WPITF  (6,2002)  I,  MI) 


P PALPHLCTION  pr  "IVAP-  AMO  "TIM"  Ic  ^•“Pr‘?S AP  Y 

n 

IF  (FLOAT  <I'(API  .Gr.0..ANn,TMAX.ST.1.  I GO  TO  40 
IF  (T-AX.PT.O.)  GO  TO  41 
THflv  = AMAxi (T1,T3>  *CT 
IF  (IVAR.GT.O.)  GO  TO  4C 
41  IVA9  = AMHV1  (OMI,OM  ll*P.*TMAXM. 

40  Of  LT=TM4X/FLOAT  ( I V/AP-  1 » 

WP1TF  (G*  20  0 21  IV!R,OrLT 
OFLT?=. 4*AMTM1( 1. /FCP.l. /PCF) 

T =OF  LT 

Y (1) = GMOL  X ( n » , 0 « I 
GO  TO  ( i PC  , 200,  3C  J)  ."'IS 

0 

c POMPPCITI-;  rILT-RS 

r 

400  ^ 410  1=1, N®’ 

410  SS(I1=S(I> 
ipf=  ic.cn 
R FAD (",FIL2» 

FF(i)=«r  B9F(?l  = ow  ?SF(MsN»  JSF(4)=RF  fSP<5)=FCr 

GO  TO  c 

4F0  J tNP  2* 2* FF ( 3 ) 

IF(J.GT.40) OO  TO  roC 
K =MP  2*  1 
pp  4 fjr  I = K,  I 
NSP=NSPf 1 

4 f-0  S(II=SS(MSP) 

IP(FF(?) .LT.RWI PW=FF(?) 

MPrJ/? 

NP’=  J 

WPITF(F>»,C01)  ( J, 9 ( J)  , 1=1, NO?) 

GO  TO  39 


C CALPMLATTON  OF  O'JTPtIT  SIGNAL  FOP  STUARE  WAVE  IN*JT  SIGNAL 


IOC  PC  101  I = 2 * T VAR 
L -UNO  <Y,  1 1 
Y (L)  =PNPLX{P.,f  .) 

TT=SMGL (Tl 

PO  1 C2  1=1, Nr»» 

X1=nFYP(-F(.J)  *TT) 

0 opjf.IT  *,MX1=",X1,-  T = ** , T T 

X?=S  (J)  ♦OMPraj 
X3=p  (J)  -OMTGAJ 
IFtTT.GT.TlIGO  TO  103 


36 


j 


non  t-  n o n 


C 


10  3 
102 


X4=0MFGAJ*TT 
«*RIMT  X4=M,X4 

Y1L)  =Y(LMAU)/?.*(<CEXP<X41-K1)  /X*MCEX«»<-X4)  -*1)/X3> 

GO  TO  10? 

YCL)=Y(L)«-Arj)/2.*X1*UCFX3  (*2*1  1)  - 1 . ) /X2  ♦ ( 0EXP(  X 3*T  1 ) - 1 . ) / X 3 > 
CONTINUE 

I F( • NOT  . I Y)  GO  T"  105 
I*1=I-1 

IF  (I. EG.?)  WRITE  (6,  2003)  IMl,Y(Htl 
WPITF  (6,2015»  I,Y<U 


CALCULATION  OF  ENIVELOP 


5 CALL  ARRAY 

IF(T.r,T,TMAX)GO  TO  EOT 
101  T=T«-OElT 
GO  TO 


CALCULATION  Op  OUT  DIJ  T SIGNAL  p 0»  COMPOSITE  WAVE  INPUT  SIGNAL 


20C  0=1. /T1 

E = 1 • / ( T 3-T  ?> 

00  211  1=1, MO? 

«m=s<n«-o*EGaj 
C (I>  =r< I) -OMEGA  I 
211  WRITC  l 6 » 30  0 4>  1,3(1), I, PfT) 

WRITE  (*=>,3003)  0,S 
OP  201  I = 2 , 1 V AR 

TT  = SNf.L  (Tj 

L =MI Nf (3, II 
Y<L»  =rMPLX<f,.,r.) 
no  20?  J=1,mo-> 

Xi=PEXP(-f <J)*TT) 

X?=CEYP(0HEGA.1*TTI 

X7=CrXP(-0MrGAJ*TT) 

X*  = l ,/C( 

X4=l./BtJ)**2 
TF(TT.Lc.T1)GO  TO  20  3 
Xf =CEXP  (O  (.1)  * TT) 
xr=CExp(r{.n*rn 
xei=r.Fxo(P(j)*Ti) 
x7i=crxo(p(j>*Ti) 

IF(TT.Lc.T?)GO  TO  204 
xf?=CFX°(p<j)*T2) 
v-'?=nr XP(P(./»  •T,» 

IFCTT.L".T3>GO  TO  2C5 
X63=CPXn  IF  ( M *T  T) 

X-’Tspr  x-»(r  { J)*T3> 

GP  TO  2<!6 

20  3 Y(L)  =Y(L»  +0*  A(J»  /'».*(  (XI- *2 1 *X4*  X?*TT<3(JI+  (X1-X3)  *X5*X3*TT/p  (J)  ) 
GO  TO  21? 


37 


THIS  NO  IS  BIST  QUALITY  PRACTICABLI 
COPY  flWHISiflBD  TODDC 


r 


i 


204  f (L>  =t(L)+n*AtJ)/-».*Xl*(IX61*(Tl/P(  J)-X4»*X71*(T1/C(JI-X5)>  f 

1 X4*X5  4-T1*  ( (X6-X61)/P(J»MX7-Xr  l>/5  C J » ) > 

GO  TO  20  2 

205  Y(L)  =Y(L>  4A(J)f2.*Xl*(0*(X44X?K61*  <T1/P(  J)-X4)4X71*  (T1/C(J)-X5) 

1 4-Tl*((XF-X61)/B(J)MX7-X71)  'C<  J>)  > -F*  (X6*  (T<B(J)-X4)  4-  X7*  IT/ 

2 C(  jy-tsi-vf}”  (T?/R(  J1-X4)  -Y^2*  (T2/C( J)-X5) -T2* ( (XS-X62)/ 

3 R(J>MXY-X*2»/C<J>n> 

GO  TO  202 

2 06  Y (L)  =y(L)4-A(J)/?.*X1*  (0*  (X4*X5*X61*  (Tl/P(  I) -X4)  f X71*  CT i/C ( J)  -X5> 

1 4-T1*  ( (XF-X61)  /B  *J)MX7-X7  1)  /C(  J))  ) 4-EMX62MT2XB  ( J)  -X41+X72* 

2 <T?fC (J)-X5)*T2*( (X6-X62) /B (J) MX7-* 72) /C  ( J)  ) -X6 3* (T 3/B ( J) - 

3 Yt,»-X77*(T3JCC  J»-X5»-T3*(  <X6-<63)  /R(  J)  KX7-X73)  /C(J)  ))  ) 

202  CONTINU* 

IF(.NOT.IY*  0,0  TO  215 
I“1=I-1 

IF  (I.F0.2)  WRITr  ce,?003)  I*l,f{111> 

WRIT  c (6,2015)  I,Y(L) 

r. 

n CALCULATION  0^  ENFV^LOP 

C 

215  CALL  ARRAY 

TF(T.GT.TMAY)no  TO  531 
201  T =T  4-  0CL  T 
o.O  TO  500 

r 

C CALCULATION  O^  OUTPUT  SIGNAL  FOP  SINF  SOJAREO  NAVE  IN=>UT  SIO,NAL 


300  no  305  1=1, NO? 
o (I) = F(I) ♦OMEGAJ 
r (i>  =s  ( I)  -one  pc  j 

305  WPITE  (6,7004)  1,3(1) , I, C(I) 

OO  301  1=2, I VAR 
L=NIN0 (3,1) 

Y(L)=C*«OLX(0.,r.) 

tt=gnol ( T ) 

DO  302  .1=1, MP? 

IF(TT.GE.Tl)GO  to  303 

Y1=C£XP(-OMEGAJ*TT» 

X2=CEXP(OWEGAJ*TT» 

R 1=T 1/°I 
X'T=0(J)*R1 
X 4=r (J) *R1 
01=° I*  TT/Tl 
n2=SIN(0i) 

03=COS(Ol) 

0t.  = X3*n2*O2-2, *0  2*03 
X *=CEXP  (-F  ( 1)  *T T ) 

X5  = ?,*PI/(B(.))*TH 

X6=X2/ (2,*XT*Y3*«,)* (04*X6) 

OL=X4*n?*r'2-2.*n2*D1 

Y7=x  1/ ( P,*X4*Xti  ♦•«.)*  044-9,  *£>!/ (2  ( J)  *T1)  ) 
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_ . 


i'A 


i * 


I 


I 


r*G£  IS  BEST  QUALlrY  praciicakle 
ISWM  CO^Y  FURMlsiiED  ioDDC 


i 


y 8«^=  x^*°T**  t ?xg=Tt**7 
R?=4.*PI»PI*T1 

Y(L)=Y(L)*A(J)MX6  + X7-X8fl/(B(  I)**?*XP«-R2*8(J))-X38/(X9*C(J)*M 
i <-r?*c(J)i) 

GO  TP 

303  X 1 = A (J) *CFXP(-S ( J ) *T  T ) *°I 

X2=«<J)*T1 

v^o  (j)  **? 

R1=T1**1 
R?=PI*PT 
P3  = R 17R2 

yt»=CPXP  f X2)  / (X3*R  *f4.*X?) 

XE=R?/(  X3*R1*4.*R?*X?> 

X2=fMJ>  *T1 
X7=C(J)**x 

X6=CFXO(X?) /(X3*R3+4.*x?) 

Y7rD2/(Y3*Rl»-4.*R?*X2> 
y (L>=Y(L)  *X  1MX4*-X6-X5-X7> 

30?  rONTINLr 

IF(.NOT.IY)  r,n  TO  3C6 
I«i=T-l 

IF  (I. £0,2)  MRITP  (6, ROOT)  IMl,r(TMl) 

WPITE  (6 , ?0 15)  I , Y (L  ) 

C - 

C CALCULATION  OF  FNFV'-LOP 

C - 

306  CALL  ARRAY 

IF(T.r-T.TM&x)r,n  TP  500 
T01  T=T*-  PF.L  T 

5 CO  I F( .MOT. IE)  GO  TO  6C0 
MPITf(6f?pn4) 

no  F0T  1=1, M 

503  WRIT r (6,2005)  I , AH ( I ) , I , T H ( I ) 

C 

C HE  TERMINATION  OF  ABSOLUTE  MAXIMUM  AMP.ITUOE  OF  ARRAY 

C 

600  CALL  MAXAMP 

WRITE  (6,2005)  AM«T 

C - — — 

C CALCULATION  OF  APRAY  CPOPDINAT«-F  FDR  CALCOMP  PLOTTER 

C - 

CALL  COORD 


C OFFLINF  CALCOMP  PLOTTING  ROUTINF 
C 

602  CALL  GRA°H(XIS,SCr,Tl,T?,T3,IGF,FF,RF) 
WRIT  F (6, VALMF) 

IF(IENO)  GO  TO  4 
2 WRITE  (6,2013) 

CALL  PLOTF(NAME) 

WRITE  (6,2014)  NAMF 


39 


o o o 


THIS  PAGE  IS  BEST  QUALITY  FRACIICABU 
ISO*  CO if  X.  EUtyUSHED  TO  DD.G  ■ ^ 


TSTOP^RFpONO  (OP) 

0FLTAT=TST03-TSTA,T 

HPTTE  TSTART ,TSTna ,n  -LT AT 

STOP 

4 IVAR-=C 
T«AY=C . 

TCP-! 

tctpp=s- nn^n (pp> 

PrLT  AT  = TSTPP-T5:T  AWT 

W PI  T r (*,1.303)  TSTAPT,TSTn-'»  ,P"LTfT 
r.r  TO  1 

7 f 0 HPlrrl& ,i*roi» 

PTOP 

4 1 3 1 FORMAT  (*0TH-;  NUMn?p  qp  polES  Of  THE  FI LTE°  EXCEEDS  THE  ALL  3 9 A TEH 
lDTMFNSION*/) 

20CC  FODM AT  (1H1,*THP  F0LL0WIN9  PORTION  0r  THIO  PRINT  OUT  IS  THE  OUTPUT 
A COTTON*) 

20C2  F PPM  AT  (1HV7//*  TVAR  = *,T7,inx,*  OPLT  = *,10D17.1J/> 

20  " 3 P ORM AT  ( 1 HI  * *TH"  C0L  LOWING  OUTPUT  IS  THE  COMPLEX  Y POR  EACH  DELT  A 
AT  TMr  t»'ORFMFNT*/7/*  = *,lP?clR.n) 

20C4  PPO*AT  ( 1 H 1 , * TH r FOLLOWING  OUTPUT  IS  THF  ARRAY  Oc  LOCAL  MAXIMUMS  ( 
ASUPSET  OF  COMPLEX  VARIABLES)  AND  TI ME  * / ) 

?0jr  CORMAT  <1H  ,*AH(*,I4,*>  - * , 1°E1  9,13  , * TH(*,I4,*)  = *,E19.10> 

290f  P CRH AT  (1H1,*THP  FOLLOWING  OUTPJT  IS  THF  ARRAY  OF  LOCAL  MINTMUMS  ( 

ACmcspT  OP  C0H°LEX  VARIABLES)  AND  TIH[*/) 

2007  FORMAT  (1H  ,*AL(*,I4,*)  = * , l^l  9. 1 G , * T»(*,I4,*)  = *.E19.10) 

20 J A POOMAT  ( 1 HC  / 7/ * ABSOLUTE  hayihUH  AMPLITUDE  IS  EQ'JAL  TO  *,1^17.10) 

2013  fcpmat  (1H0////*  ’'HI  S STATtm^JT  TERMINATES  THIS  RUN*) 

2P14  format  (1H0//7*  THf  NUM3FR  OF  n.OCKS  ON  THE  3LOTTINC  ta»E  IS  *,IF) 
23  lc  FORMAT  < 1 H ,*Y(*,I6,*)  = *,10?E19.1C) 

33„"  FORMAT  ( 1HC ,*P<* , 12,*)  =*, 1P2T19. 13) 

3001  FORMAT  (1M0,*S{*,T2,*)  =* , 1 p?£1 9 • 10  » * S(*,12,*)  =* , 2E1R. lu> 

3002  FORMAT  (1H0,*A(*,T2,*)  =* , 1R2E19 . 13 ) 

3003  FORMAT  (1HC,*0  - * * 1R“19. 10  , * i =*,E1R.10) 

3004  FOPM/ir  (1H3  ,*9(*,T2,*)  =* , 1P2C1 9 . 13  , * 0(*,I2,*>  = *,2EH.10) 

4000  FORMAT  (1HC,*THIS  RUN  STARTFQ  AT  TIME  *,F7.3,*  AND  “NO ED  AT  *,F7.3 

A,*  TAKInc,  * , ft,  3 , * SECONDS  TO  EX  ECUTf*  ///) 
cup 

surpoutinf  array 


T . I • i 1 4, * ) = **E19,10) 


S ( * , 12  * * > =*,2E19.1u) 


ROUTINF  TO  LTMTT  0'JT°UT  OF  FNVFLOP  ARRAYS 


complex  Y<3) 
nrUPLF  PELT , T 

0 TMF  NS  I ON  AH(?000l ,TH(?orO) 

COMMON/VAPS/AH,  AMPT,Q5'LTtTVAR,L*M,Tf  TH,TMAX,XWAR,XVAR1, 
AYVAR  *Y ,OELTt , ifr  , TFA 
IF  < L# E o»  2)  CO  TO  ? 

Ir  (L.NF.3)  CO  to  7 

YMAX  = AMAX1  (PEAL  (V(L-P))  ,RF  AL  (Y(.-l)  ) t REAL  (Y  ('.)  )) 
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IF  <TMAX.EO.RF4L (»<L-1> > » GO  TO  3 
r,p  to  F 

3 IFfY (L-’l .EO.rMPLXfO. ,0. )» GO  TO  f 

M=H*1 

IF(M,r,T.  2C0r ) GO  TO  6 
TW(M|=SN0l (T) 

ah<m)=real<t<l-i>> 

IF(M,lE.  ■»!»  GO  TO  F 
T (L-U  =y  (l>  rOMPLY  <0., 0.) 

T=T»0ELT2-0ELT 

C HPITE<6,3000»Y(L-1I,TH<M) 

GO  TO  f 
2 M = 1 

TH(M)=SNGL (T-OELT) 

AM<mrRLAL<Y<L-l)» 

RFTURM 

5 T <L-2»  =Y  (L-1T 
Y<L-D=Y(L> 

PFTUR*' 

6 WRITE  (G,2000> 

R TOP 

7 WRTT c <S,2Q01» 

STOP 

2000  FOPMAT  <1H1,*THC  "LOTTING  A PR  AY  rXGffOR  THE  At  L00  ATEO  Dimension  rv 
A 200  0*1 

2001  FORMAT  (1H1,*FRP0"  IN  T'<f  C AL 1UL  A’’!  ON  nr  L COUNTER*) 

3000  FORMAT<fx,*y=*#io-»e1R.10,yJ  ^A<  r = *,E10.13» 

3001  FppMAT  <FX,*Y  = *f  1P’E1R.10»*J  MIN  T**,f  1A.MI 

FMO 

SUBROUTINE  MAXA-" 

C 

C ROUTINE  TO  CALCULATE  A3SOLUTc  MAX  SIGNAL  A MOL I TUOE 

C - 

COWPLrX  Y<3> 

OOUBLC  OrL T ♦ T 

DIMENSION  AH(2C0C)»TH(2J00) 

COMMON/VARS7AH,AMc»T,nFLT,IVAP,L,  m,t,TM,TMAX,*VAR,WAR1, 
AVVAR,Y.OFlT2,IF", IF A 
AMPT  = ah ( 1 ) 

no  l i="*m 

1 AMPT  = AMAX1  ( AM®T,  AM(in 
RrTURM 
FT'O 

EIJOROUTINF  COGRO 


C ROUTINE  TO  CALCULATE  PLOTTER  0 OOROI  MAT " S I”  TNO-rs 

C 

COMPLrX  Y<T» 

DOUBLf  OrL  T . T 

nlMFNEIOM  AM(?000) ,TH(?0C0l 

COMMON/ VARS /AM, AM"T ,DCLT , TV AP  ,L,M,T,TH,TMAX,<V4?,XVARl, 
AT VAR  »T,0cLT,,IFo,TFA 


■ r--  v 
• » - 


I 
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YVAR1=X'»AR 

2 IF(«:-Lr'AT(*'-n.wc..^VeRl*10C.)G'i  ro  1 

XVAP1=XMAR1*1.1 
GO  TO  ? 

1 HPI T E (F..20C0)  XVAR1,YVAP 

2300  FORMAT  (lH0///<*  LENGTH  OF  THE  X-AXIS  IS  * , c . Z , * INCHES,  AN  0 THr 
AWIGht  of  th?  Y-A*IS  IS  *tFfi,i,*  IYCw^S*) 
fmo 

fuorO'JTTMF  GRAPH  (yIS,SCF,T  1 ,T?,T  3,IOF, FF, SF) 


P0UTTNF  to  PLOT  ARRAYS 


C0MPLrY  y<») 
nnu°Lc  0cl  T , T 
LOOTOAL  IFO,IFA 
0 IMF  fjF  I O*1  FF(F»  , Sr(f  > 

OIMENSION  AH(  ?G00>  ,TH  (230?) 
niMFNSIOM  IT(A)  ,ITL(5) 

COMM ON/ VARS /AM,  AM«T,OcLT,IVAR,L,  M,T,  TH,TMAX  ,xvar,xvari  , 

1Y  VAR ,Y , OEL  T?,IF°, TFA 

0 AT  A (IT  (I),  1 = 1, 8)  <6MSOtlARc,  OHOOioosi  TE,  lnHSINE  S0UAR,2HE0, 

1 1CHRUTTFRN0RT  , 1HH,  RHCHF^r  SME‘1,  fhpolES/ 

OAT A (I TL(I)  ,1  = 1,51/6,9,12, 11, 1/ 

AH(Mf i) =TH (MfH =0. 

= TMAX/XVAP1 
AH{H«-?»  =AMoT/?.S 

call  AX  IS (3 . ,-YVAR,15HTIMt  (MSEC  0^0 S ) , -15 , X « AR1 , 3 . , TH(M+l),TH(Mf?> 
1) 

T YV= YVAR*?, 

rALL  AXIS  < 3. ,-YVAR, 18HPFLA TIVR  L I T HO?  , 1 A , T Y\f , RO  . , - 1 . , . 4 > 

CALL  PLOT  ( 1 . , C . , 31 
CALL  PLOT  (XVAR1, 

CALL  °LOT  (o.,p,,T| 

CALL  LINF(TH,AH  ,m,i,o.  ,«*> 

00  2 1 = 1, M 
2 A H (T  > =- AM ( I ) 

CPU  LIN?  (TH,ah,m#i,  o,,4» 

XPLHS=XVAP14-3.?5 
C TmpuT  SIGNAL 
h = .1  4 

C ALL  fymool  ( J.E,- 3,7,H,  lXMINP’IT  SIGNAL  t , 0. , 13) 

CALL  Symbol  -3.7  ,H,  IT  ( XT  3) , 0 . , I TL  (^IS)  ) 

OALL  SYMPOL  (2.37,-3.R,H,  3hf  ^,0.,3I 

CALL  fymroL  (?.!»/+, -3.  9,  .ct,  iHC,C.  ,1) 

FC=S0C 

CALL  HUMBFR(2.77,-3.0,H,Fn,  0.,2> 
o PL  L SYMBOL  M.ET,-3. g, h,thmh7,0.  , ?» 

CALL  SYMBOL ( 2.3T,-4. 1 ,h, tmT  1=,3.  ,t> 

YT=T1 

CALL  NUMBER  (2.77  ,-4. 1»H,  XT , ".,3» 

CALL  SYMBOL (3. 5 7, -4. 1 , H , 4HM SE 3, C . ,4) 
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If(T2.F,'.C  ,)r,0  Tfl  5 

PALL  PY'inOl  (2.37,-4.3  ,H,  3HT  2=,0.  , 3> 

YT  = T? 

PALL  ►•lJHr»FR  (2.7-'  ,-4.  3 ,Hf  YT»C  . ,3) 

CALL  ^YHOCL  <7.5  7,-4. 3,H,4HMSrC,0  . ,4> 

I F (T  3. E°  • C . >00  Tn  5 

CALL  SYMBOL  (2.37,-4.5,H,3HT3=,0.  ,3) 

*T=T  7 

CALL  MU^BFR('».7-',-4.5,h,yt,  C.,3) 

CALL  SYM°OL  <3.F7t-it,£>,H,4HMScS,3.»4> 

C FILTER 

k CALL  SYMqCL(4.E,-'».7,H,7uriLT'£«t  , C.  ,7» 

K^F^tl  > 

r'PLL  SYMBOL  (S,6,-7.7,M,IT  (?*K’c«-3)  » ] • * I T L f ^r+J)  ) 
PALL  SYMBOL  (5.6 ,-7.9,H,7Mr  =,].,3) 

CALL  f'YMPOL(l3.7  3,-3.  9 ,.  07,lHP,0.  ,1) 

CALL  NU^«PR(1>.i,-3.9,H,Fr(F)  , 3.» 

PALL  5YMPOL(7.1»-3.o,H,  3HMH7,').»  3) 

PALL  NUM0FR(5.*,-h.1,H,FF(3>  , D.,-l» 

CALL  SY-POL (5.9,-fe.l,H,IT{«)  , ).,5) 

CALL  ^vmool  (5.6t-4.3,H,3HPW=,  1.,  3) 

CALL  NI»MOFR(<j.l,-4.3,H,FP(?)  , 2) 

PALL  PYMOpL(7.i,_u.3,Ht3MMH7,"!.»  3) 

IP(FF(4)  .F0.3.»p,0  TO  A 

CALL  symbols. f>, -4. ?,M,5HRF  = ,C.,5> 

PALL  .R(6.?,-4.5,H,FF(4»  , 1.,  “>> 

8 T F ( I PF , E 0 , i ) CO  to  1G 

Xr7.?  »YY=-7.7 

CALL  py-TOL  <X,XY,M,7MPILTTR»  , 7) 

VrYfl.i  *YP=S*TCU 

CALL  SYMnPL(X,YY,^tIT(2*KFf 3) ,Q.,IT_<<F«-3)> 
vy=y Y-. ? 

PALL  SVMflPL TX,XY,Hf3HF  =,P.,7> 

PALL  Symbol  (X  «•.  1 3,XY,.C7,1HC,Q.,1) 

fp=SF(5) 

CALL  NUMPf  R(Xf.«5,XY,H,PO,0.  ,7* 

CALL  PYMOOL  (X4-1.  E »XY  ,H,  7HMHZ,  , 3) 

XVsYY- . 7 

CALL  NIIMHE  R (X,XY  SF(3I  , A.  ,-i) 

CALL  SYMBOL  (X*-.  7,  YY,H,IT  (8)  ,C.,5) 

XY=XY-. 7 

CALL  c‘YMBOL  (X,XY  ,m,  7HBM=t  0.  , 3> 
fp= Sc(2) 

PALL  NUmBFR(X4-.F,XY,H,FC,J.  ,2> 

PALL  SY1«OL<XH.P,XY,M,3H*1H7fJ.,  3> 

IF(SF(M  ,FP.  0.)  CO  TO  10 
y y=y y-«  ? 

CALL  PY1POL(X,XY,^,3HRFr,0. ,3» 

CALL  MllMRFRfXf.*,tY,H,SF<4)  ,C.,2) 

1C  IF(.NOT.  IFO»r,0  TO  40 
I M=m_e 
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40 

* * 


AH  (Mf  !)=••, AM(Hfl) 


00  20  1=1, I* 

AH(1 1 = 1,  / CTH  ( IfH  — TH<  11  >-SCF 
AH(W)  = AM<H«-il=AH<  10) 

1 HsM-l 

no  30  1=10, IM 
AH(M)=AHIN1 (4M(M) ,AH(I)I 
AH<m*-1)^ AHAY1IAM(*1h-H  ,AH<I)  > 

AM(Mfl)  = (AM(HH)-AH(M))/9,5 

ORINT  AM(H)  =-,ah<H) 

TM(M)rp. 

TH(MM)=TMAX^XVAR1 

CALL  AXT^(XVARl,f>,,15MFRtn'J«:N-;Y  IM  H Z, -15 , Y V AR  , 9 0 . , AH  ( M>  , A>Lt  HU)  > 
CALL  DASHLN(TM(10)  ,AHU0)  , I M-q,l) 

TF(.N0T.IFA)G0  TO  40 

WRITE (6, 50  CI,AM(I)  ,I,TH(I)  ,T  = 10  ,IH) 

FORMAT  (*1  THT  FOLLOWING  OUTPUT  IS  30M'»UTr0  SIGNAL  FRFQ UE^CT * 

1 (*  FP(*ti4t#>  = *,E19.10,*  TH<*,I4,*)  = *,F.19.10>> 

CALL  °LOT<X°LUS,0«t-3) 

RFTURM 

END 


I 


